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ABSTRACT 

In a disk around DM Tau, previous observation of 13 C0 ( J = 2 — 1 and 1 — 
transitions) derived the 13 CO gas temperature of ~ 13 — 20K, which is lower than 
the sublimation temperature of CO (20 K). We argue that the existence of such 
cold CO can be explained by a vertical mixing of disk material. As the gas is 
transported from a warm layer to a cold layer, CO is depleted onto dust grains 
with a timescale of ~ 10 3 yr. Because of the steep temperature gradient in the 
vertical direction, an observable amount of CO is still in the gas phase when the 
fluid parcel reaches the layer of ~ 13 K. Apparent temperature of CO decreases 
as the maximum grain size increases from //m-size to mm-size. 

Subject headings: stars: planetary systems: protoplanetary disks — stars: pre- 
main-sequence — ISM: molecules 



1. Introduction 



Theoretical models previously showed that the protoplanetary disks can be divided to 



three layers in the vertical direction from a chemical point of view (lAikawa fc Herbst 111999 



Willacy fc Langer Il2000l ; lAikawa et al. 1120021 ; iBergin et al. II2006I ). Surface layer of the disk 
is a PDR (photon-dominated region) because of the interstellar UV, stellar UV, and stellar 
X-rays. As we go deeper into the disk, the radiation becomes weaker, and the density gets 
high enough to allow high abundances of organic molecules. In the midplane, on the other 
hand, heavy-element species are mostly depleted onto grains because of very high density 
and very low temperature. The three-layer model naturally explains the characteristics of 
the disk chemistry observed in radio wavelengths: low average abundances of gaseous organic 
species and relatively high abundances of radical species (ex. CN/HCN ratio) in comparison 
with the molecular clouds JDutrev et al. Ifl997l : loi et al. Iliooa iThi et al. 1 120041 ) . Molecular 
emission lines arise mostly from the intermediate molecular layer, because it is chemically rich 
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and spans a density range that covers critical densities of many (sub-) millimeter rotational 
transitions. 



Dartois et al. I (120031 ) observed two rotational transitions (J — 2 — 1 and 1 — 0) of various 
CO isotopes, from which they derived 12 CO temperature to be 30 K and 13 CO temperature 
13 — 20 K. Since the 13 CO lines trace the deeper layer of the disk than the 12 CO lines, it 
was the first observational evidence that the gas temperature increases with height from the 
midplane at radii of r > 50 AU, which are traced by the current radio telescopes. In addition, 
the 13 CO temperature challenges the theoretical model; since the sublimation temperature 
of CO is ~ 20 K, a layer with T ~ 13 — 20 K is considered to be the freeze-out layer. 



Recently, lAikawa fc Nomura I (120061 ) briefly argued that existence of such cold CO can 
be explained by vertical mixing due to turbulence. In the vertical direction gas can migrate 
over ~ 10 AU in the vertical direction, and thus can travel from the T = 20 K layer to the 
T = 13 K layer within 10 4 yr, which is comparable to a freeze-out timescale in the DM Tau 
disk. In this paper, we examine their arguments in more detail using various disk models. 
The rest of the paper is organized as follows. In §2, we derive a typical migration length of 
material within the freeze-out timescale of CO. In §3, a distance between the T = 13 K and 
20 K layers is measured in various disk models, and is compared with the migration length. 
We compare our results with the recent numerical calculation of turbulent disk chemistry in 
§4. Summary is contained in §5. 



2. Basic Equations 



Following I Willacy et al. I (120061 ). we consider a vertical mixing in the disk. The net 
transport flux of species i by the turbulent diffusion is 



-Kn{E 2 



dx{i) 
dz 



where x(i) is the abundance of species i relative to the number density of hydrogen molecule 
n(i) / n(H.2) ■ We assume that the diffusion coefficient K is equal to the viscosity coefficient: 
K =< v t l >= ac s h where v t is the turbulent ve locity, I is the mixing length, a is the alpha- 
viscosity parameter (IShakura fc Sunayev Ill973l ). c s is the sound velocity, and h is the scale 
height of the disk. The chemical continuity equation is 



dn(i) d(j)(i) 



Pi 



dt dz 

where Pi and Li are the chemical production and loss terms for species i. 



(2) 
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We now consider the steady-state abundance distribution of CO near the boundary 
between the warm molecular layer and freeze-out layer. Since CO is the most abundant 
and stable species among the carbon-bearing species in the warm molecular layer, the main 
production and loss terms are desorption from and adsorption onto grains, respectively. 
As the matter is transported from the warm molecular layer to the freeze-out layer, the 
desorption term becomes negligible. Hence the vertical distribution of CO is determined by 
the diffusion from the warm molecular layer and the adsorption onto grains, i.e. 



d^{CO) 2 8kT 

= -na Si 77^ n g n ( CO ) ( 3 ) 

dz y 7rm(CO) s 

where a is the grain radius, S is the sticking probability, k is the Boltzman constant, m(CO) 
is the mass of CO molecule, and n g is the grain number density. In order to derive an 
analytic solution, we neglect the density gradient (dn(H-2) / dz) , which is justified because the 
abundance gradient is steeper than the density gradient. Then the basic equation is reduced 
to 

d 2 x{CO) 2c / 8kT 

K — j-^ = na 2 S x — — n g x CO). 4 

dz 2 y 7rm(CO) s 

The analytic solution for this equation is 

x = x exp ( 5 ) 

where x is the canonical CO abundance (10~ 4 ) in the warm molecular layer and z is the 
height of T = 20 K layer, i.e. the boundary between the warm molecular layer and freeze-out 
layer. The decay length A is given by 

x fna 2 S I 8kT \ ^ 

A= ^V^(co) ng J • (6) 

In other words, material can migrate from the boundary z = z to z = z — A within the 
freeze-out timescale of CO. 



3. Results 

In this section, we evaluate A in various disk models. Considering the vertical temper- 
ature gradient, we discuss if the turbulence can transport observable amount of CO to the 
T < 20 K layer. 
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3.1. Migration length of CO and vertical temperature gradient 



The decay length A is determined by the ratio between the rate coefficient of adsorption 
and diffusion coefficient. The former depends on temperature, sticking probability, and total 
grain surface area. Temperature is set to be 17 K, since we are interested only in a region 
between 13 K and 20 K, and since the dependence of A on temperature is not strong. Sticking 
probability S is set to be 0.5. As for dust grains, we assume dust/gas mass ratio of 1% and 
grain radius of 10~ 5 cm. The number density of grain n g is then 1.6 x 10~ 12 nH, where nn 
is the number density of hy drogen nuclei. If the grain size distribution follows the MRN 



model (IMathis et al. 1 119971 ). the total surface area (J 7ia 2 n g (a)da) of dust grains should 



be about three times larger than assumed in our model. Howeve r, observations suggest 



that the dust grains in disks are larger than the interstellar dust (jD'Alessio et al. 1 12001 



Rodmann et al. II2006I ). Dependence of our results on grain size in disks will be discussed in 
the next subsection. 

The diffusion coefficient depends on a and scale height. We a s sume a = 0.01. Analyzing 
the CO emission lines in DM Tau disk, iGuilloteau fc Dutrey I (119981 ) found the intrinsic 
local velocity dispersion i s essentially thermal, with a turbulent component of 0.05 km s _1 
< v t < 0.11 km s -1 , while iDartois et al. I (120031 ) estimate v t ~ 0.16 km s -1 . Thus by setting 
a = 0.01, we assume the mixing length I of several 10~ 2 h. The scale height of the disk is 
given by 



h 



2kTr s 



(7) 



GM*[im H 

where /i is the mean molecular weight, mn is the hydrogen mass, G is the gravitational 
constant, and M* is the mass of the central star. We assume /x = 2.34 and M* = O.5M . 
Then A is given as a function of hydrogen number density Wh and disk radius r: 



A = ( 1.6 x 10" 



12 



2na 2 S / GM*/jm H 
ac s y 7rm(CO)r 3 



-1/2 



n H 



Here we evaluate A in two disk models from literature: D M Tau disk model, wh ich is 
constructed to fit the SEP assumi ng the grain size of 0.1 /im (ICeccarelli et al. 1 120051 ) . and 
a model by lD'Alessio et al. I (119991 ). who assumed interstellar-type dust is well mixed with 
gas. Models give nn and T as a function of r and z. At each radius we calculate A using 
% in the layer of T = 13 — 20 K. In Figure 1, asterisks, open boxes, and crosses depict A 
with the density rin at T = 20 K layer, T = 13 K layer, and with the average density of 
T = 13 — 20 K layer at each radi us, respectively. The b lack marks are the DM Tau disk 
model and the gray marks are the iD'Alessio et al. I (119991 ) model. For example, % is ~ 10 7 
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cm 3 at the boundary between the molecular layer and the freeze-out layer at r ~ 400 AU 
in both models. Then the decay length A is ~ 10 AU. 

Also shown in Figure 1 (closed circles) is the distance between the layers of T = 20 K 
and 13 K, 6^20-13, at each radius in the two disk models. It can be seen that at radius of 
> 400 AU, (IZ20-13 is about 4 times larger than A. In other words, when the disk material 
transported vertically from the 20 K layer to the 13 K layer, gaseous CO is e -4 times the 
canonical value. 

It should be noted that the column density of CO gas with T < 20 K is given by 
An (CO), where n (CO) is the number density of CO molecule at the 20 K layer (z = zq). 
At the radius of 400 AU, for example, A is about 10 AU, and no (CO) is ~ 10 3 cm -3 . Then 
the column density of 12 CO with T < 20 K is ~ 10 cm , which corresponds to the opacity 
of r ~ 1 for the J = 1 — emission line of 13 CO, assuming the isotope abundance ratio 
of n( 12 CO)/n( 13 CO)= 60. Therefore, the amount of col d 13 CO produced by the turbulent 



mixing is high enough to account for the observation by iDartois et al. I (120031 ) . 



Finally, it is obvious from the equation (6) that in the disks with larger (smaller) a, A 
is larger (smaller), and thus the cold CO is more (less) abundant. 



3.2. Effect of grain growth 

So far, we have assumed that dust grains in disks are similar to the interstellar grains. 
However, observations suggest grains in disks are larger than the interstellar grains. Dust lane 
width of edge-on disks is smaller than that of model disks with interstellar grains. The SED 
i s better reproduced by the disk models with the maximum grain radius in millimeter range 



(ID'Alessio et al. Il200ll ). In this subsection, we compare the decay leng th A and the vertical 



temperature gradient in the disk models of lAikawa fc Nomura I (120061 ) . They investigated 



the effect of grain growth on disk structure (% and T); the number density of grains with 
radius a is assumed to be proportional to a -3 5 , and the maximum size a max is varied with 
the dust /gas mass ratio fixed. 

Figure 2 compares A and 0^20-13 i n the models with (a) a max = 10/zm and (b) 1 mm. 
Compared with the case of uniform grain size (0.1/im), the total surface area of grains is 
smaller by a factor of (a) 0.3 and (b) 0.03. Because of the smaller total surface area of 
grains, the adsorption time scale and thus A are larger in the model with a max = 1 mm (see 
eq. 6). On the other hand, ^20-13 is smaller in a max = 1 mm model than in a max = lO/xm 
model at r > 100 AU. Temperature in the disk atmosphere is lower, and hence the vertical 
width of the disk is smaller in the model with larger a max , because of the reduced number of 
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small grains which mainly contribute to the absorption of stellar light. At r = 100 AU, the 
temperature in the midplane layer increases with grain growth, and dz2o-i3 is slightly larger 
in the a max = 1 mm model. We can conclude that at r > 100 AU the amount of cold CO 
gas increases as the a max increases from 10/xm to 1 mm. 



Discussion 



Recently ISemenov et al. I (120061 ) investigated the molecular evolution in a turbulent disk 
by coupling the chemical reaction network and turbulent diffusion. They showed that the 
two-dimensional diffusion increases the amount of cold CO, but the amount is too small 



to account for the observation by iDartois et al. I (120031 ) . Why is the cold CO much less 
abundant in their model than estimated in the present work? 

We suspect that the deficiency of cold CO is caused by the too efficient grain-surface 
reactions. Even in the T ~ 20 K region in their static model, CO abundance is much 
smaller than the canonical value. Inclusion of the vertical diffusion reduces the cold CO, 
because more CO is converted to H 2 CO by grain-surface reactions (see their §3). Since the 
sublimation temperature of U2CO is higher than that of CO, the hydrogenation depletes CO 
at T ~ 20 K. The rate of hydrogenation, however, seems to be overestimated in their model, 
as discussed below. 

In the layer of T ~ 20 K at r > 100 AU, we estimate the number density of H atoms 
(not included in H2 ) in t he gas phase to be n(H) < 10 cm -3 , referring to the model of 
Aikawa fc Nomura I (120061 ) with a max = 10/im. The number of H atoms accreting onto a 
grain with radius a is 



8kT 



n(HW = 2 x 10- 4 ( — ^ 
v ; V20K7 



1/2 



7rm(H) 

The desorption rate, on the other hand, is 



0.1/zm 



-1 



(9) 



2n s ff D 
7r 2 m(H) 



exp 



Ed 
kf 



(10) 



where n s is the number d ensity of adsor ption sites on a gr ain (~ 1.5 x 10 15 cm" 2 ), an d Ed 
is the adsorption energy (IHerbst I Il993l ). If Ed is 350 K (IHasegawa fc Herbst Ill993l ). the 
desorption rate is 7.4 x 10 4 s _1 at 20 K. Hence a number of H atoms on a dust grain is much 
smaller than unity even for grains with a = 10/im. With such a low H atom abundance 
on a grain, the hydrogenation rate is limited by the accretion rate rather than determined 
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by a migration rate on a grain (iCaseUi et al. Ill998l ). According to the model description 
of ISemenov et al. I (120061 ). they seems to have estimated the hydrogenation rate using the 
migration rate, which can be an overestimate. 

Reality can be more complex than the above argument. For example, grains are not 
spher ical but would be poro us, which will enhance the possibility of reaction before desorp- 
tion (IPerets fc Biham 1 120061 ) . On the other hand, laborator y experiment suggests t hat CO 
is hydrogenated only in the surface layer of the ice mantle (IWatanabe et al. 1 120031 ) . Most 
of the currently-available astrochemistry models may overestimate the hydrogenation rates, 
since they do not discriminate the surface species from species inside ice mantle. Modeling 
the grain-surface reactions is still a challenging subject in the astrochemical studies. 

However, it is an observational fact that cold (~ 13 K) CO gas exists i n the disk. If 
the CO depletion is caused by hydrogenation in the model by lSemenov et al. I (120061 ). it tells 
us even the abundance of simple gaseous species can be dependent on how we model the 
grain-surface reactions - an important caution for future models. 



5. Summary 

We have shown that in protoplanetary disks the turbulent diffusion can produce observ- 
able amount of CO that is colder than the sublimation temperature. As disk material is trans- 
ported to the layer with T < 20 K, freeze out proceeds with timescale of ~ 10 3 (nH/10 7 cm~ 3 ). 
Since the gas migrates over a distance of ~ 10 AU within this timescale, CO gas of column 
de nsity ~ 10 17 cm" 2 exi sts in the cold (T < 20 K) layers, which explains the observation 



by iDartois et al. I (120031 ). As the dust grains in the disk become larger by coagulation, 
the freeze-out time scale increases, and the vertical temperature gradient at T ~ 20 K gets 
steeper. Hence, the apparent CO temperature decreases as the maximum grain size increases 
from /zm-size to mm-size. 
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Fig. 1. — The black marks depict DM Tau disk model by Ceccarelli. et al. (2005), while 
the gray marks are the disk model by D'Alessio et al. (1999) with the accretion rate of 
1O~ 8 M yr _1 . Asterisks, open boxes, and crosses depict A with the density n H at T = 20 K 
layer, T = 13 K layer, and with the average density of T = 13 — 20 K layer at each radius, 
respectively. In the DM Tau disk model, for example, the densities wh at T = 13 K and 20 
K layers are 6.6 x 10 6 cm -3 and 1.2 x 10 6 cm" 3 at r = 616 AU, while they are 4.2 x 10 8 cm -3 
and 9.1 x 10 7 cm -3 at r = 191 AU, respectively. Closed-circles depict the distance between 
the layers of T = 20 K and 13 K {dzza-is)- 
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Fig. 2. — Decay length of CO abundance A and distance between the layers of T = 20 K 
and 13 K (^20-13) as i n Figure 1. Disk models are by Aikawa & Nomura (2006) with the 
maximum grain size of 10 /zm (a) and 1 mm (b). In the model with a max = 10/im (a), the 
densities wh at T = 13 K and 20 K layers are 2.4 x 10 7 cm" 3 and 1.2 x 10 7 cm" 3 at r = 305 
AU, while they are 6.8 x 10 8 cm -3 and 3.7 x 10 8 cm -3 at r = 100 AU, respectively. In the 
model with a max = 1 mm (b), the density n-^ at T = 13 K and 20 K layers are 1.9 x 10 7 
cm -3 and 1.4 x 10 7 cm" 3 at r = 305 AU, while they are 5.2 x 10 9 cm -3 and 2.9 x 10 9 cm -3 
at r = 50 AU, respectively. 



